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Abstract 

We derive and study quasicanonical Gibbs distribution function 
which is characterized by the thermostat with finite number of parti- 
cles (quasithermostat). We show that this naturally leads to Tsallis 
nonextensive statistics and thermodynamics, with Tsallis parameter 
q is found to be related to the number of particles in the quasither- 
mostat. We show that the chi-square distribution of fluctuating tem- 
perature used recently by Beck can be partially understood in terms 
of normal random momenta of particles in the quasithermostat. Also, 
we discuss on the importance of the time scale hierarchy and fluc- 
tuating probability distribution functions in understanding of Tsallis 
distribution, within the framework of kinetics of dilute gas and weakly 
inhomogeneous systems. 

1 Introduction 

Nonextensive statistical mechanics introduced by Tsallis [|l[] has been devel- 
oped and applied by many authors H]-[B|; see M for a recent review. Particu- 
larly, various aspects of the nonextensive statistics and thermodynamics have 
been studied by Abe and Rajagopal [Q. Derivation of power law canonical 
distributions from first principles of statistical mechanics has been analyzed 
by Almeida ||; see also |L0 . 
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The entropic parameter q of Tsallis statistics is associated to the measure 
of nonextensivity, with the usual extensivity recovered at q = 1. In Tsallis 
formalism, q appears in the definition of entropy which generalizes the usual 
Boltzmann-Gibbs-Shannon entropy. 

In this paper we show that a nonextensivity characterized by a parameter 
naturally arises from the consideration of the microcanonical Gibbs distribu- 
tion. Namely, in Sec. ^| we consider a subsystem weakly interacting with a 
system containing big but finite number M of particles which we call qua- 
sithermostat. It appears that the resulting equilibrium distribution function 
and the associated entropy (Sec. |^) reproduce exactly those of Tsallis statis- 
tics, with q — 1 ~ 1/M. In the (thermodynamic) limit, M — > oo, the usual 
canonical Gibbs distribution is recovered. 

Recently, Beck || have studied application of Tsallis formalism to turbu- 
lent flows and achieved a very good agreement with experimental measure- 
ments. Particularly, it is worthwhile to emphasize that no fitting parameters 
have been used by Beck to reproduce histogram of the acceleration of a test 
particle advected by the turbulent flow, while some other authors have used 
four free parameters in some ad hoc modified exponential function, for this 
purpose (see [|] for details). The main idea underlying Tsallis statistical me- 
chanics approach to turbulence is to introduce fluctuation of temperature or 
fluctuation of energy dissipation [[7], || described by x 2 distribution of order 
m. 

In this paper we show that the relevance of x 2 distribution can be un- 
derstood in terms of normal random momenta of particles in the quasither- 
mostat, with m = 3M representing the total number of random variables 
(Sec. |). 

Some aspects of fluctuations and time scale hierarchy in dilute gas and 
weak inhomogeneous systems are presented in Sec. []. 

In this paper we show the importance of the time scale hierarchy and 
fluctuating probability distribution functions in understanding Tsallis distri- 
bution, which is viewed as the average value of fluctuating canonical Gibbs 
distribution over the random parameter f3. 

2 Quasicanonical Gibbs distribution 

We start with the classical microcanonical Gibbs distribution, fi(x,a, E), 
which describes isolated equilibrium system of L classical particles with total 
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energy E, 

h(x, a, E) = —±—5(H(x, a) — E), (1) 

where f2(a, E) = J 5(H(x, a) — E)dx, H(x, a) is Hamiltonian of the system, 
a denotes external parameters, dx = IIk^l^, and Xi = (pi, r*) are phase- 
space coordinates of the particles. 

In the closed system of L particles we select a subsystem of N 3> 1 
particles which is surrounded by the subsystem of L — N = M particles. As- 
suming M 3> N, neglecting the interaction energy between the subsystems, 
and taking the limit M — > oo, from the microcanonical distribution ([J) one 
can derive, following one of the known procedures, the usual canonical Gibbs 
distribution, 

Hn(x, a 



f N (x,a,T) = -^exp 



kT 



(2) 



where Z = J exp [—Hn(x, a)/kT] dx is the classical statistical integral. Here, 
k is Boltzmann constant and T is an absolute temperature of the subsystem 
of M = oo particles (called thermostat); dx = lli<i<jv^ x i- 

It is worthwhile to note that the canonical Gibbs distribution comes also 
as the only physically relevant solution of supersymmetry condition in the 
context of topological quantum field theory approach to Hamiltonian sys- 
tems |T]. 



Below, we follow some procedure assuming that the number of particles 
in thermostat is finite, M < oo. Such a thermostat can be referred to as a 
quasithermostat. 

The total Hamiltonian is represented by 

H(x, y) = H N {x) + H M {y) + H- mt (x, y), (3) 

where x and y denote phase-space coordinates of particles in the considered 
subsystem and in the quasithermostat, respectively, and H- mt (x,y) character- 
izes interaction between the subsystem and quasithermostat. Herebelow, we 
drop dependence on the external parameters to simplify notation. 

Since the total system is closed, it is described by the microcanonical 
distribution, 

f N+M (x,y,E) = C8(H{x,y) - E), (4) 

provided that it is normalized, J fN+Ai(x,y, E)dxdy = 1; dy = rii<i<Af dyi- 
Note that the numbers M and N are additive integrals of motion. To obtain 



3 



the distribution function describing the considered subsystem one should 
integrate out phase-space coordinates y of the quasithermostat, 

f N (x) = CJ5(E- H N (x) - H M (y) - H- mt (x, y))dy, (5) 

For big N and weak short-range interactions one can ignore the interaction 
energy, H int (x, y) <C Hn(x), in the above expression. By representing Hamil- 
tonian of the quasithermostat as a sum of the kinetic and potential energy, 

H M (y) = K M (p 1 , . . . ,p M ) + U M (n, . . . ,f M ), (6) 

where we take for simplicity 

M -2 

and Um{t) = Sj^-$(|ri — fj\) for point-like particles of mass m, we can 
rewrite (||) as 

f N (x) = cj6(E - H N (x) - U M (r) - z 2 )5(z - K\{ 2 )dzdpdr, (8) 

where new variable z has been used; dp = Ili<i<M dpi, dr = IIi<i<M dFi. 
Performing the integral over z with the use of 

j8(z-K)l\p))dp^z* M -\ (9) 

and general property of Dirac delta function, 

/ 5{g{z))dz = f £ 5i f~ Z ^ dz, g{z k ) = 0, (10) 

we get, up to normalization constant, 

f N (x) =CJ(E- H N (x) - UMir))^- 1 dr. (11) 
We rearrange this as follows 



3M -i 
2 ± 
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3M , 

x (E-U M (r))-- l dr } (12) 

where we have denoted, for brevity, H{x) = H N (x). 
Denoting the kinetic energy per particle by 

3 _ K M (p) 
2(3 ~ M ' 

which tends to constant for big M, we obtain that for M ^> N, 

E - U M {r) ~ — , (14) 

which does not depend on r. 

Inserting (|14|) into ([12]) and performing the integral over r we obtain, up 
to normalization constant, 

3M 1 

M*) = C (l - m PH{x)) 2 • (15) 

We note that to obtain the canonical Gibbs distribution (Q) it remains to 
take the (thermodynamic) limit, 



3M 1 

Jfo^ (l - -j3ff(x)j 2 = exp [-/3if(z)] (16) 



and identify = fcT. 

However, we are interested specifically in the case when M is finite so 
that we do not take the above (thermodynamic) limit, and turn to analysis 
of the distribution function (p"5|). 

For M > 1 we make the approximation 3M/2 — 1 ~ 3M/2. Then, 
denoting 

<?-l = ^7>0, (17) 
H 3M ~ ' y ' 



we have from ( 15Q 



f N (x, a) = ±- (1 - (q - l)PH{x, a))T* 



= ^-^ q _ x [-pH{x,a)}. (18) 
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Here, we have introduced Z q = J exp q [—j3H(x, a)]dx due to the normalization 
condition, / fN<ix = 1, g-exponential is defined by expjz] = (1 + sz) 1 ^, 
s = q — 1, and we have restored dependence on the external parameters a. 
In view of the fact that the number M is taken to be big but finite, i.e. 
q ~ 1, the normalized distribution function flT8| ) can be referred to as the 
quasicanonical Gibbs distribution. 

The obtained distribution function (|IBD strongly resembles the distribu- 
tion function in Tsallis nonextensive statistics 0, |2|, where q is called 
entropic parameter. 

We remind that g-exponential, exp„ Jz], reduces to the ordinary expo- 
nential, exp[z], at q — > 1, so that the canonical Gibbs distribution is recov- 
ered. 

Typical plots of exp^-J— z] at various q are presented in Fig. |l]. In the 
physical context, to avoid the negative and increasing values of exp (? _ 1 [— z\ 
one must restrict positive z by a suitable condition (the high-energy cut-off). 

Some useful formulas with g-exponential are: 

dexf> q _ x [±z\ldz = ±(ex Pq _ 1 i±z}) 2 -«, (19) 
J exp^zjdz = ±q- 1 (ex Pq _ 1 [±z]y, (20) 

exp q -i[-z] = (exPi- q [z]) > ( 21 ) 
exp^ [z x ] exp g _ x [z 2 ] = exp g _ x + z 2 + (q ~ l)ziz 2 ] , (22) 
the g-Poisson integral, 

z^e-lidz = z^ ai^^l 1^.^.(9- ^ (23) 



and 



o 



z n e~lfdz = he - cg)' (n+1)/2 x 



for n > —1, 9+ 2(iIg) +1 > 0) ^2^" — T^q > ^' an< ^ rea ^ c ' wnere 2-fi is hyperge- 
ometric function and T is Euler gamma function. 

The reflection of the g-axis with respect to the point q = 1, maps exp (J _ 1 [z] 
to exp-^_ q \z\ — (1 -\- (1 — 9)2) which is often used in the literature as a 

definition of the g-exponential. 
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Figure 1: The g-exponential exp -J— z] = (1 — (q — l)z) 1 ^ q x > at various 
values of q. 

Also, it is remarkable to note that the slightly generalized ^-exponential, 

ex Pg _ lib [z] = (1 + (q - l)z)^ + \ (25) 

can be represented as a power of the ^-exponential, 

(1 + ( ? _ = (exp 9 _ l [z)) e , (26) 

for the parameter c = 1 + (q — 1)6. For example, (1 + (q — = 
(l + (g— l)z) q ^ q ~ 1 ^ = (exp Jz]) 5 , so that for small g — 1 we have exp J,?] ~ 
(exp q _ 1 [z]) q . 

Note that (exp.-Jz]) corresponds to the scaling of z in the q — > 1 limit, 
namely, (exp g _ 1 [^]) c — > exp[cz], while in general (exp Jz]) ^ exp^Jcz]. 

In contrast to the canonical case, the distribution function ( p~8|) depends 
not only on the temperature (/3 _1 ) but also on the total number of particles in 
quasithermostat displayed by the parameter q. Also, the property (|22| ) of q- 
exponential leads to a non-factorizable character of the distribution function 

Using the identity exp -J^z] = exp{lnexp -Jz]} and the new variable 
z q = lnexp (? _ 1 [^] = (g — l) _1 ln(l + (g — l)z), one can formally rewrite g- 
exponential in the usual exponential form, exp^-Jz] = exp[z g ]. This allows 
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to reformulate any g-exponential distribution function, e.g., flTB|), in a fac- 
torizable form with respect to z q . Consequently, one can formally achieve 
an additive character of Tsallis formalism with respect to such "g-deformed" 
entities z q . 

Almeida |§ has derived Tsallis distribution (|18|) from the first principles 
of statistical mechanics using the parametrization 

d9 . 

5I = «- 1 ' < 27) 

where 9 = 1/(3 is the statistical temperature. Tsallis parameter q is given a 
physical interpretation in terms of heat capacity of the heat bath due to the 
fact that 9 is related to the temperature of the heat bath. For an infinite 
heat capacity (thermostat), q = 1, one recovers canonical Gibbs distribution 
while for finite heat capacity (quasithermostat), q 1, one gets power law 
distributions. 

The micro canonical approach meets this viewpoint as it establishes de- 
pendence of q on the number of particles M of the quasithermostat which 
evidently characterizes its heat capacity, d9/dE ~ 1/M. It is worthwhile to 
note however that in general the form of the dependence 9(E) in the l.h.s. of 
Eq. ( ^7|) seems to be not fixed by first principles of statistical mechanics. The 
linear dependence, 9(E) = (q — 1)E, does yield Eq. (^7|) while a quadratic 
dependence would require yet another parameter. In contrast, the used mi- 
crocanonical approach with the kinetic energy in the form (0) suggests only 
one parameter, M, which we relate to q via Eq. (|T7|). 



3 The entropy 

The entropy associated to the distribution function f N = f N (H(x,a), (3,q) 
defined by ( |T8"D can be readily obtained, up to a factor and constant term. 
First, we identically rewrite Eq. (|18|) as 

f N = exp q -i\PF q -(f3-(q- l)E q )H(x,a)}, (28) 

where we have denoted 

F^^ViV' (29) 
and g-logarithm is an inversion of g-exponential, hx s [z] = (z s — l)/s, s = q—1. 
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Some useful formulas with g-logarithm are: 

limlng.x z = \nz, ln ? _i z~ x = — ln^g z, (30) 

ln ? _!Z p =pln p( g_i)Z, z ln^-i z = z q hii_ ? z, (31) 
dln 9 _i = z^ 2 = + (g - 1) ln ? _! z), (32) 

y ln g _i zdz = g~ 1 z(ln g _ 1 z — 1), (33) 

and 

ln 9 _i zxz 2 = ln 9 _i zi + ln g _i z 2 + (g - 1) ln g _i Zl ln 9 _i z 2 , (34) 
Then, with the use of the inverse function, 

H = ptf N ) = (PF q - ln ? _! f N )/(p - (g - (35) 

the normalization condition, / f^dx = 1, the thermodynamic heat, 5Q = 
J H(x,a)(5fisr)dx, and the integrability condition, dQ = T(f3)dS q one can 
obtain (cf. 0) 

S q (a, T) = -J f N m g _i f N dx = -J f q N lni_ g /Web. (36) 
Here, the representation 

rfN 

S q = T (0) / (p(z) + const) dzdx (37) 
Jo 

with a suitable constant has been used. The integrability factor is found to 
beT(f3) = (q[3-q(q-l)F q )-\ 
With f3~ x = kT, we obtain 

6Q = {w ~ q{q " 1)Fq ) dSq - (38) 

In the discrete case characterized by microstates with energies E n , the 
obtained quasicanonical distribution function (|18D and the entropy (|36|) are 
written as 

/ n = ^-exp,_ 1 [-^ n ], (39) 



n=l n=l 



Z q = J2exp q _ 1 [~(3E n } J £/ n = l, (40) 
9 



w w 

S q =-J2fn H-l fn = -T, fn ^l- q fn (41) 
n=l n=l 



V W ( fl — f \ X^ W fl _ 1 

_ ^n=lUn J"/ _ Z^n=l Jn x 

g — 1 q — 1 

The entropy ([IT]) coincides with Tsallis entropy [|TJ, with g treated as 
an entropic parameter, and reduces to Boltzmann-Gibbs-Shannon entropy, 
S = — fn hi f n , at q — ► 1. It is remarkable to note that either /„ and 
can be used for the averaging in Eq. (pEID provided a suitable choice of q. 

It is well known that by extremizing the (convex) entropy (|41"| ) under 
suitable constraints || one obtains Tsallis distribution (|39|). 

Due to the property ( |34f) of g-logarithm the entropy (^l|) is non-additive, 

5) = S 9 (A) + S q (B\A) + (q- l)S q (A)S q (B\A), (42) 

where we use the conditional entropy for statistically dependent systems A 
and B (e.g., subsystems with long-range interactions) based on the condi- 
tional distribution function f(A, B) = f(A\B)f(B) = f(B\A)f(A). 

In the case A and B are statistically independent in a conventional sense 
(e.g., subsystems with short-range interactions), we have formally S q (A\B) = 
S q (A) and S q (B\A) = S q (B). 

We note that the (q — 1) term in Eq. ( |42"D is non-zero even for such 
"statistically independent" systems contrary to a naive expectation that the 
statistical independence should imply an additive character of the entropy 
S q . The origin of this implication seems to be due to the fact that here the 
parameter q does not depend by construction on a scale or correlation param- 
eter, such as a correlation length, which allows one to determine a statistical 
independence of the systems A and B in a physical context (when the sys- 
tems A and B are separated by distances much bigger than the correlation 
length their evolutions are physically independent). 

Formally, for a constant parameter q ^ 1 the systems appear to be al- 
ways statistically dependent on each other. Deeper analysis is required for a 
consistent formulation of the statistical independence within the framework 
of non-additive formalism ([42]). 
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4 Normal random momenta and fluctuating 
statistical temperature 

Recently, Beck H have studied application of Tsallis formalism to turbulent 
flows and indicated that the factor of Tsallis equilibrium distribution function 
( |39| ) can be obtained by averaging Boltzmann factor of the canonical Gibbs 
distribution ||, 

J e-fc-fiftdp = (1 - (q - l)/3E n )T*, (43) 
o 

where 1 — q = 2/m and 

m ~ 

1 ( n \ 2 ~ m n md\ . A A . 

(44) 



' m ' 
2 



m/3 



2/5 



is the probability density of the x 2 distribution of order m (the distribution 
of sum of squares of m normal random variables with mean zero and unit 
variance), with the average (p) = J °° $f(/3)d/3 = [3 and the relative variance 
((/? 2 ) — (f3) 2 )/(P) 2 = 2/m. This is viewed as the result of fluctuation of 
temperature or fluctuation of energy dissipation ]7|, ^ and has been found of 
particular relevance for nonequilibrium systems. Beck have used an assump- 
tion that the time scale on which (3 fluctuates is much bigger than that of 
kinetic relaxation of the system. This provides a quasistatic behavior of the 
system with respect to variation of j3. 

We see that x 2 randomization and averaging of the canonical Gibbs dis- 
tribution over random f3 leads to Tsallis distribution. 

It seems that the relevance of \ 2 distribution ( JH| ) to relate Boltzmann 
and Tsallis factors is not due to a coincidence. By comparing this approach 
0, H with the above obtained results flTTD and flTB|), we see that the num- 
ber of random variables is m = 3M, which is exactly the total number 
of components of three-dimensional momenta Pi of M particles constitut- 
ing the quasithermostat. Since the quasithermostat is assumed to be in a 
quasi-equilibrium state we can treat the momenta pi (microscopic character- 
istics) to be normal random variables with mean zero and unit deviation as 
a simplest model. Hence, the kinetic energy (^|) is x 2 distributed due to its 
definition. We associate the kinetic energy of particles in quasithermostat to 
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its statistical temperature in accord to Eq. flT3|). Hence, becomes x 2 
random variable. 

This partially clarifies the physical origin of the x 2 distribution used by 
Beck in Eq. (|43|). However, we note that in our approach \ 2 distribution 
arises for the statistical temperature (3~ l (not for the inverse statistical tem- 
perature (3). 

Usually, normal (not x 2 ) distribution of fluctuating temperature with the 
relative variance from Eq. (|58|) is used as a good approximation. We note that 
the averaging of the canonical Gibbs distribution over normally distributed T 
(i.e., (3~ l ) with the mean T and some deviation a is not analytically tractable, 
while the averaging over normally distributed (3 yields the analytic result, 

fe'^giPW = e~^ 2H2 (l - erf(^^)) , (45) 

where 

^) = ^=ex P [-%^] (46) 

is normal distribution, and erf(x) is the error function. The resulting expres- 
sion (fHI) includes one special function and is less handful as compared with 
Tsallis distribution fl4~3|). 

Since the entropy S q and mean energy H of the equilibrium system are 
not additive for q 1, the function F q is not additive as well. Indeed, from 
Eqs. (|T8|) and ([36|), we have 

S q = (i3-(q-l)F q )H-i3F q , (47) 

where H = J Hf^dx. This implies that (3 is not a universal parameter of 
the system, in contrast to the case of additive entropy and mean energy. 
Indeed, f3 is not the only parameter characterizing the quasithermostat, with 
q being the second temperature parameter measuring the relative variance 
of the fluctuating parameter around the value (3 viewed as an average. 

From the viewpoint of canonical Gibbs distribution, accounting for a 
finiteness of thermostat made in Sec. automatically arises to the averaged 
X 2 distributed fluctuation of the parameter (3 around the mean (3 with the 
relative variance ~ 1/M. Actually we associate the fluctuations of temper- 
ature with eventual fluctuations of the energy of a finite thermostat. We 
conclude that the x 2 statistics is implicit in Tsallis statistical mechanics. 
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The spatial scale dependence of the entropic parameter, q = q(r), intro- 
duced by Beck |5j can be readily understood since q is related to the number 
of particles. For bigger spatial scale the fluctuations become smaller so that 
q —>■ 1 in the big scale asymptotics. 

To illustrate how the resulting nonadditivity affects the temperature char- 
acteristics, we consider two systems with Hamiltonians Hi and H 2 , respec- 
tively, at the same (3. Abe and Rajagopal |12| have elaborated the power-low 
distribution based thermodynamics of equilibrium to much extent. Below, 
we discuss this issue in the context of fluctuations of the parameter (3. 

For weakly interacting canonical systems we have exp[— j3H\\ exp[— /3H 2 ] = 
exp[—/3H], where H = Hi + H 2 is conserved. As usual, this means that 
the two equilibrium systems compose an equilibrium system at the same 
(3, while at different /3 there is no equilibrium, with conserved H. The 
X 2 randomization, (3 — > (3, and averaging of each of the two canonical 
systems, with the same average j3 and relative variance q — 1, gives us 
exp -J— exp q _ 1 [—/3H 2 ] = exp^J— [3H e g], where we have denoted H efi = 
H+ (q — l)f3HiH 2 . On the other hand, the averaging of the composite system 
yields a different result, (exp[— /3(Hi + H 2 )])p = exp q _i[— (3H]. 

It is misleading to interpret the above relations as that the composite 
system appears to be not at thermal equilibrium despite the identical (3. A 
thermal equilibrium of two systems in the sense of their absolute tempera- 
tures is well formulated for a non-fluctuating (3 but not for the average, /3, of 
fluctuating (3. The extra /3-term in H e g looks like a pure effect of the x 2 aver- 
aging since it is proportional to the relative variance q — 1. Indeed, in general 
{f0,H x ))p{f0 L H 2 ))p ± {f0,H x )f0,H 2 ))p regardless to factorizability of 
the function f(f3, H) with respect to additive H. Even more difficulties in 
understanding thermal equilibrium would arise when one naively considers 
two systems with the same average (3 and different relative variances, q\ — 1 
and q 2 — 1. 

In order to clarify the physical content of the problem of fluctuating (3 we 
turn to the time scale hierarchy emerging from the microscopic consideration 
of a dilute gas and weakly inhomogeneous systems of a finite characteristic 
size. 
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5 The time scale hierarchy 



Internal characteristics of a dilute gas with density n = N/V are known to be 
the interaction range r int , the mean distance between particles, r avc ~ n -1 / 3 , 
and the free-flight length r fli ~ l/(nr i 2 nt ). The main parameter of a dilute gas 
is the density parameter, 

e = nrf nt < 1. (48) 

For example, ambient air is characterized by the typical value e ~ 10~ 4 . For 
small e we have scale hierarchy, 

n n t < r avc < r fli , (49) 

with rf nt : r 3 vc : r|; = 1 : e^ 1 : e~ 3 , from which it follows the time scale 
hierarchy, 

Tint < 7- avc < Tffl, (50) 

where r int = r int /t> is the interaction time, r ave = r ave /v is the average time, 
Tfli = Tfli/f is the free-flight time, and v is some characteristic velocity of 
particles. Typically, for dilute gases r int — 1CT 12 sec and r n ; — 1CT 8 sec. 

For such systems the quasi-closeness and quasi-equilibrium of subsys- 
tems are well defined, and a description in terms of one- and two-particle 
probability distribution functions, or mean local densities (particle number, 
momentum, and energy), is rigorously justified. 

One-particle distribution function varies much slower than the higher or- 
der ones due to the smallness of e, providing thus quasi-conservation of the 
mean local densities. At the microscopic level a quasi- additivity of the (quasi- 
) local energy density (hence, of the energy) can be seen from the last term 
in its definition, 

m = E *)R72m + E *(\* ~ rM (51) 

i j^i 

Here the last sum makes the biggest contribution for the coordinates fj of 
particles which are at distances less or about r- nit from the point r. In fact, 
whether h(r) can be taken additive or non-additive depends on the scale. 
We start by considering a sequence of systems, 

Ni C N 2 C • • • C Ni C • • • C Woo, (52) 

generated by the relations between the numbers of particles in these systems, 
1 < Ni < N 2 < • • • < Ni < • • • < A^oo = oo. The systems M x and can 
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be viewed as a smallest system and an ideal thermostat, respectively, in the 
equilibrium state. More specifically, the systems M% are treated as thermody- 
namical systems of characteristic sizes which are external characteristics, 
such that 

Tint < r avc < r x < ■ ■ • < r fc _i < r fli < • ■ 

■ • ■ <^ r k <ti ■ ■ ■ <ti Too = oo. (53) 

Here, we have used the fact that for a dilute gas the free-flight length is much 
bigger than the interaction range, ^> r int , so that some number (k — 1) of 
systems fall into the range between r ave and r fli . In this paper, we will not 
consider such systems which require a separate study 

Actually we assume that Moo can be represented as a collection of weakly 
interacting equilibrium systems Mi so that all bigger systems M%, % > 1, are 
in incomplete equilibrium state at some time scale. 

We remind that in the framework of the canonical Gibbs distribution, the 
interaction between different Mis is taken negligible, £ = H-^/H^ = (zero 
order approximation in small £). However, the relaxation of some system Mi, 
i > 1, to an equilibrium state is solely related to H int so that at least first 
order approximation in small £ should be used to describe the relaxation 
process. 

The kinetic relaxation length r re i ~ r^, and the associated kinetic relax- 
ation time, T re i ~ rfli, are much bigger than the correlation length, r cor ~ r int 
and the associated correlation time, r cor ~ r int , for e <^ 1. 

In terms of small parameter e, one can define (see e.g. [|13|| ) 

ri = y/erm, n = y/erm, Ni = l/\/e, (54) 

for the smallest system Mi- The size t\ = vt\ is taken to provide conser- 
vation of the local densities of particle number, momenta and energy in the 
approximation of weak inhomogeneity, 

Tcor/n < 1, (55) 

at the time scale t\. For e = 10~ 4 , we obtain the minimal number of particles 
Ni — 100. To specify the hierarchy, we can put Ni = 1/e'^ 2 . 

To each characteristic size r^, i > k, we can associate the gas dynamical 
relaxation time, 

n = ri(ri/r fli ) 2 , i > 1, (56) 
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so that the time scale hierarchy is 

Tint < Tavc < T fli < T k < T k+1 ■ ■ ■ < Too = OO. (57) 

For % > k, the sizes of the systems Mi are much bigger than r ni so that 
for such systems one more small parameter can be introduced, Ki = rm/ri, 
Knudtsen parameter, in terms of which n = T\/Kf. 

It should be emphasized that due to N\ = 1/ \fe the e = approximation 
in the (nonffuctuating) canonical Gibbs distribution implicitly means that 
Ni = oo (the approximation of a continuous medium for Ni particle system). 
Evidently, the consideration of bigger systems Mi, % > 1, becomes of no sense 
in the e = approximation. This is the point where a large scale dependence 
is totally ignored within the approximation of nonfluctuating distribution 
function. To keep this dependence one should derive the quasi-equilibrium 
distribution function in at least the first order approximation in e and allow 
distribution function to fluctuate. 

6 Fluctuating probability distribution func- 
tions 

Now we turn to fluctuations. Following to Onsager principle, kinetics of a 
linear non-equilibrium thermodynamical system (its relaxation to an equilib- 
rium state) defines fluctuations in the system. In turn, the fluctuations of the 
thermodynamical entities (evolution of the time correlation functions of their 
deviations from equilibrium values) can be equivalently described in terms 
of Langevin equation with 5-correlated stochastic sources. Thus the fluctua- 
tions of equilibrium systems can be viewed as a stationary process with the 
time scales defined by the relaxation times r», i > k, of the gas dynamical 
stage of the relaxation. As shown in Sec. |5] these time scales are known to 
be much bigger than that of the kinetic relaxation r re i thus providing a local 
equilibrium framework. 

We are interested in a finite thermostat so that we take some finite size 
equilibrium system A/} serving as a quasithermostat for the system Mi, i.e., we 
assume Nj ^> N^. In accord to the above consideration, the time scale Tj of 
the quasithermostat Mj at which its f3 (kinetic energy per particle) fluctuates 
is much bigger than the time scale Tj at which any thermodynamical entity 
of the system M% fluctuates. In general, there is finite number of big systems 
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■M+i, J\fi + 2,---,J\fj, each making contribution to fluctuation of the tempera- 
ture with different relative variances at corresponding time scales r i+1 ,...,Tj. 
Observed marginal distribution for M% will depend on the observation time 
T bs as compared to these time scales. We remark that the observed relative 
variance of fluctuating (3 is thus depends on the scale and r b s - 

At the time scale r b s ~ t% the canonical Gibbs distribution makes a 
very good description of A/i, with (3 treated as a constant, and an additive 
character of the entropy and mean energy. 

At long time scale, r b s > Tj ^> r*, the fluctuations of (3 can not be 
ignored so that Tsallis-like distribution emerges as an effective description, 
i.e., we should replace (3 — > (3, and make the averaging {e~ ,3Hi )p 1 with (3 = {(3) 
and some relative variance. This results in a non-additive character of the 
entropy and mean energy of A/i- Note however that in the latter case the 
usual quasi-equilibrium (thermal equilibrium) still holds for any given value 
of j3 at the small time scale r,. This picture meets that assumed by Beck ||. 

For j — > oo (an ideal thermostat), we have Tj —>■ oo so that j3 is constant 
at any time scale r > 7*, thus q — 1 for A/^. 

To put the above consideration on solid grounds, below we outline statis- 
tical mechanics approach to fluctuations, within a dynamical framework. 

We can treat the quasithermostat A/} as an equilibrium system, which 
can be described in pure additive thermodynamical terms, with the canonical 
Gibbs distribution, Z~ x exp[— Hj/kT], providing the averaging in J\fj. 

However, kinetic and thermodynamic fluctuations naturally occur in the 
quasi-equilibrium system A/} and can be calculated in a linear approxima- 
tion (small fluctuations) for the case of weak inhomogeneity of the system, 
Eq. fl55|). We notice that rigorous description of fluctuations is out of scope 
of the canonical Gibbs distribution and the one-particle distribution function 
entering Boltzmann equation. 

In a dynamical framework, one should consider Bogolyubov equations 
with a short-range interaction potential, $(|r m — f n \), for point-like particles 
(providing convergence of the integral for local energy density), small den- 
sity parameter e (i.e., three- and higher order particle collisions are ignored), 
and weak inhomogeneity (providing small gradient of the one-particle distri- 
bution function), and account for the two-particle correlation function with 
a partial relaxation of initial correlations. This implies fluctuations of the 
one-particle distribution function; see e.g. and review fH| . The form 
of the resulting kinetic equations allows one to study separately implications 
related to Brownian effect (a discrete collision structure of the medium) thus 
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discarding influence of fluctuations on kinetic processes. 

The resulting picture is equivalent to consideration of some Langevin 
equations with 5-correlated stochastic sources, from which one can extract 
two-point correlation functions and spectral densities of the gas density, ve- 
locity, and pressure, for a local equilibrium state. In the local equilibrium 
approximation correlation functions of the stochastic sources can be found 
also by phenomenological methods of non-equilibrium thermodynamics |15 



The spectral densities can be used to estimate the gas dynamical re- 
laxation time, 1/r ~ vk 2 ~ nk 2 , where v is kinematic viscosity, k is heat 
conductivity, and k is wave vector. For a finite system of characteristic size 
L, the wave vector k is naturally bounded from below, \k\ > 1/L, from which 
with the use of estimation v ~ k ~ vr^ it follows Eq. (p^). 

We focus on the temperature fluctuations, for which one can obtain the 
approximate variance, 

((AT) 2 ) ~ kT 2 (dT/dE) v , (58) 

with temperature T being normally distributed (in this approximation), and 
T in the r.h.s. is given in terms of equilibrium value. This means that 
T fluctuates about the average T with the relative variance k/C v , where 
C v = (dE/dT)y ~ Nj is a heat capacity of the quasithermostat Nj, thus 
ending up with the consideration of Sec. |4|. 

It is worthwhile to note that the entropy S of Nj fluctuates as well, 
{(AS) 2 ) ~ kC p , C P ~N. 
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Clearly, natural fluctuations occur not only in the quasithermostat Nj 
but also in the considered smaller system N%. However, we can ignore these 
fluctuations as they make small contribution to specific fluctuation of the 
temperature T governed by the quasithermostat, due to AT* <C Nj. 

We point out that a major difference from the above consideration arises 
when turning to fully developed turbulence considered in |5j . Beck used sta- 
tistical temperature related to kinetic coefficients (fluctuations of which are 
evidently not related to a finiteness of quasithermostat) and have averaged 
the probability distribution function (describing a stationary non-equilibrium 
state) derived from a nonlinear Langevin equation which is different from 
the canonical Gibbs distribution. Fluctuations of the friction coefficient (or 
viscosity) are seem to be of a kinetic type originating from large scale fluctu- 
ations of the phase-space density function averaged over the smallest scale. 
Such fluctuations may have a big relative variance as they are not related 
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to the particle number of some big system, with the main contribution to 
the energy dissipation happening at the smallest scale being fluctuations, 
Hi ~ 5ui, where u, is velocity and 5ui is fluctuation of the velocity. Indeed, 
for Kolmogorov (the smallest) length scale 77 of the turbulent flow Beck has 
derived under simple assumptions a theoretical value of the entropic param- 
eter, q = 3/2, which essentially deviates from unity (the relative variance is 
considerable, q—1 = 1/2). One can see || that Tsallis formalism does work 
in this case as well. For bigger scales, fluctuations of the energy dissipation 
decrease so that in zero order approximation in rj/r the energy dissipation 
is constant. Clearly, one should assign the value q = 1 for this large scale 
asymptotics. 

Fluctuations of f3 in the canonical Gibbs distribution mean that we deal 
with a fluctuating distribution function. The \ 2 average value of this fluctu- 
ating distribution function is Tsallis distribution, 

(e-* H b = efff- (59) 
The relative variance of this fluctuation is 



where s = 1 — q. For small (1 — q), this relative variance behaves as (1 — q) 3 . 
Also, up to a first order in (q — 1), we obtain 

exp^l-PH] ~ (l + q -^f3 2 H 2 ^} expi-PH}. (61) 

It is remarkable to note that this function is non-factorizable with respect to 
H because of the presence of H 2 . 

From the above dynamical consideration it follows that fluctuations of a 
probability distribution function arise naturally in studying both equilibrium 
and non-equilibrium phenomena. However, in studying equilibrium state one 
usually uses non-fluctuating distribution function, thus discarding specific 
effects such as wellknown Brownian one. Namely, the e = approxima- 
tion means a continuous medium approximation, at which Brownian motion 
can not be described, even for the smallest physical equilibrium system Ai 
consisting of, e.g., Ni = 100 particles for e = 10~ 4 . 

In contrast, the use of a phenomenological Tsallis approach makes pos- 
sible to account for some of the effects related to natural fluctuations of 
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an equilibrium distribution function or that corresponding to a stationary 
non-equilibrium process. 

Finally, we point out that even small fluctuations of the one-particle dis- 
tribution function, 5f/f <C 1, make an essential contribution (as big as 
Boltzmann integral which is treated as a leading term) to Boltzmann equa- 
tion with Langevin sources [[13]. This indicates importance of accounting for 
fluctuations of probability distribution functions, both for equilibrium and 
non-equilibrium states, at some scales. 



7 Discussion 

The microcanonical approach in the form presented in this paper evidently 
can not cover a universal character of Tsallis formalism. First, since we as- 
sume M ^> 1 the parameter q — 1 is very close to zero that allows us to use 
expansion in q — 1 and calculate small corrections to the canonical Gibbs 
distribution, while for considerable q — 1 the arguments of our derivation in 
Sec. § break down. Second, we have ignored the interaction energy H int that 
is not valid for systems with long-range microscopic interactions to which 
Tsallis statistics is applied. Also, the nonextensivity can arise not only as 
the effect of a finite heat bath. For example, Tsallis distribution arises in 
the description of polydispersed colloids with the fluctuating parameter be- 
ing particle size |J, and was used to describe statistical properties of fully 
developed turbulence |7|, |8j ; see also review in JJ . 

In general, fluctuations of the parameter (3 (or statistical temperature) 
and averaging of the fluctuating distribution function over x 2 -like distribution 
can be viewed as a defining feature of Tsallis-like nonextensive statistical 
distributions, regardless to the origin of fluctuations. 

Adopting this viewpoint, one can try various types of distribution of j3 
in replace of x 2 distribution, and allow considerable relative variance of /3. 
For example, it is natural to assume Poisson distribution for the random 
momenta Pi in replace of the normal distribution; see also [TO. Also, one 
may be interested in a more detailed description and account for higher 
order momenta such as biquadratic deviations of (3. 

From this point of view, Tsallis formalism can be treated as a suitable 
universal tool to deal with fluctuating distribution functions without turning 
to underlying kinetic issues. However, it should be noted that the use of 
an averaged fluctuating distribution function means incomplete statistical 
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description. 

Finally, we stress that the very notions of extensivity (additivity) and 
intensivity in thermodynamics are essentially based on the requirement that 
the system is homogeneous, which is provided for big systems with weak 
interactions or, more precisely, in the thermodynamic limit, N, V — > oo, 
N/V = const. These notions make no strict sense for inhomogeneous sys- 
tems such as finite systems or systems characterized by the size of the order 
of correlation length r cor . As a consequence, in this case there is no ther- 
modynamical equivalence between equilibrium microcanonical and canonical 
ensembles. 
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